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Abstract. - We investigate the quantum ratchet effect under the influence of weak dissipation 
which we treat within a Floquet-Markov master equation approach. A ratchet current emerges 
when all relevant symmetries are violated. Using time-reversal symmetric driving we predict a 
purely dissipation- induced quantum ratchet current. This directed quantum transport results 
from bath-induced superpositions of non-transporting Floquet states. 
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■ An intriguing phenomenon in non-equilibrium trans- 
. port is the ratchet effect [1-3], i.e., the emergence of 
' directed motion in the absence of any net bias. Net 

transport results from an interplay between ac driving, 

spatio-temporal asymmetries, and non-linearities in a pe- 

. riodic potential. This mechanism provides the basis for 

■ an increasing number of experiments ranging from particle 
transport in biological systems [4] and nano-engines [5] to 
charge transport in semiconductor heterostructures [6,7], 

. superconductors [8] and spin transport [9] . Symmetry in- 
' vestigations revealed the necessary conditions on the ac 
. force and the static potential, such that a ratchet current 
' can emerge [10-13]. 

A widely employed model for studying the ratchet ef- 
, feet is a one-dimensional periodic potential in which clas- 

■ sical Brownian particles move [1-5]. It describes also 
\ the motion of a thermal cloud of cold atoms in an ac- 

■ driven optical potential [13]. As the atom cloud is cooled 
' down further, one expects quantum effects to become rele- 
vant [14]. The Hamiltonian limit of such quantum ratchets 
has been studied recently [15-17]. A more realistic de- 
scription of quantum ratchets necessitates inclusion of the 
ubiquitous decoherence and quantum dissipation [18-21]. 
For moderate-to-strong dissipation, incoherent tunneling 
transitions prevail and the quantum ratchet current can be 
studied within quantum rate theory [18-21], while in the 
high-temperature limit, one can employ a Fokker-Planck 
equation with quantum corrections [22]. For very strong 
friction, a description in terms of an effective Smolu- 
chowski equation comprising leading-order quantum cor- 
rections is appropriate [23]. By contrast, the crossover 
towards the coherent quantum regime, i.e., the under- 



damped regime [24], in which already weak decoherence 
significantly alters the Hamiltonian dynamics, still repre- 
sents an ambitious challenge. 

In this letter we study ac-driven quantum ratchet trans- 
port in the technically demanding regime of weak quan- 
tum dissipation where quantum coherence and relaxation 
affect each other. We analyze within a Floquet-Markov 
description [25] the dynamics on quantum attractors by 
expanding them into the Floquet states of the correspond- 
ing coherent time-dependent system. Then a most intrigu- 
ing question is whether violation of time-reversal symme- 
try due to weak quantum dissipation is perceivable in the 
quantum attractor and in the quantum ratchet current. 

Model and quantum master equation. A quan- 
tum particle in a time-dependent periodic potential obeys 
the Schrodinger equation 



.. d , , , 
^ri—tp{x,t) = 






V{x,t) = Vouix)-xE{t), 

where u{x) = u{x + L) with max]u(a;)] ~ 1 describes 
the shape of the L-periodic potential. The driving E{t) 
is a time-periodic field with zero mean, E{t + T) = 
E{t), {E(t))T = 0. Henceforth, we use I/^l = L/2Tr, 
{m/kj^Vay/^, and Vq as units of distance, time, and en- 
ergy, respectively, such that formally k^ — m — Vq = 1, 
while h becomes the effective Planck constant ?ifcL/^/"^Vo 
[14]. 

By the gauge transformation \tp) -^ exp{—j-xA{t))\ip), 
we bring the Schrodinger equation (1) to the spatially pe- 
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riodic form [17] 
d 



in 



dt 



^{x,t)^(^[p-A{t)f+u{x))i,{x,t), (2) 



with the vector potential A{t) = — L E{t')dt' and the 
momentuni operator p = —ihd/dx. The corresponding 
Hamiltonian has recently been realized in cold atom exper- 
iments [14]. The Schrodinger equation (2) is time periodic 
with period T = 2tt/ll! and, thus, according to the Floquet 
theorem, it possesses a complete set of mutually orthog- 
onal solutions of the form [^^(t)) = e~"^°*/''|(/)a(t)). The 
Floquet states |(/'a(t)) ~ {(pait + T)) and the quasienergies 
Ca, —huj/2 < Sa < tiiu/2, are obtained from the eigen- 
value problem [H (t) - ihd/ dt]\(l){t)) = e\(j>{t)) [26]. Owing 
to discrete translation invariance, all Floquet states are 
characterized by a quasi momentum k with \(f>a(x + 2TT)) — 
exp{iTiK)\(j>a{x)) . We restricted our study to states with 
K ~ 0, which can be expanded into the plane waves 
\n) = (27r)~^/^exp(inx). Physically, these states corre- 
spond to initial states with atoms populating all wells of 
the spatial potential equally. Such initial conditions are 
the natural ones for the recently proposed [27] ring-shaped 
optical potentials which have already been realized exper- 
imentally [28]. 

We incorporate decoherence and dissipation by coupling 
the driven system (2) to a bath of non-interacting har- 
monic oscillators [29]. Following a standard approach to 
weak quantum dissipation [25] , we decompose the reduced 
density operator g into the Floquet basis of the coherent 
system, Qafjit) — {(\)a{i)\Q{i)\'i)(i(t)) ■ Assuming that dissi- 
pative effects are relevant only on time scales much larger 
than the driving period T, we arrive at the master equa- 
tion 



QaH = —T^^a — ^fj)gaf3 + 2_^ i^afi.a'fi' Qa'fi'i 

a'P' 



(3) 



with the time-independent transition rates 

n 
~ ^PP' /_^ NjS'i a' ^nXafj" ,-nX p" a' ,n f j\ 



P",n 



Oaa' 7 ^ Na" p' ^nX p' a" ,-nXa> 



(3,n ■ 



where X„/3.n = ((0„(t)|xe-'"'^*|(/.^(t)))r and 7V„;3,„ = 
N{ea - e/3 + nhuj) with N{e) = {■je/h^)nth{e)- Here 
'^th(e) = [exp(e/fcBT) — 1]~^ is the thermal occupation 
number and (• • • )t denotes the average over one driving 
period. 

This Markov approximation requires that the coupling 
strength 7 is the smallest frequency scale in the prob- 
lem, such that 7 <C k^T/T% and 7 <C Aap/h, where 
Aa/3 == jfa ~ e/3| IS any splitting of the Floquet spec- 
trum [25]. The latter condition is rather strict because in 



the present case, the quasienergies are even dense on the 
interval [—huj/2,fiuj/2]. Thus, the condition 7 <C A.ap/h 
is violated for any finite dissipation strength. Yet it is 
obvious that only a finite number of Floquet states signif- 
icantly contributes to the density matrix at large times. 
We thus validate our results for the asymptotic state g^p 
with the following criterion: We sort the Floquet states ac- 
cording to their weights Qaa and consider only results for 
which the first Ng states fulfill the condition j/Aap < e, 
where s is the threshold value. 

For any initial density operator gap{0), the solution of 
the master equation (3) converges to a unique "quan- 
tum attractor" being the fixed point g'^p = gap{t -^ 
00) of the quantum master equation. Note that in the 
Schrodinger picture, the density operator is periodically 
time-dependent and, thus, describes a limit cycle. Since 
in Floquet representation, the attractor gi° ^ is neverthe- 
less time- independent, the asymptotic current, defined as 
the time-averaged momentum expectation value, reads 



•^ = X! SlpPap; Pap = {{(j)a{t)\p\(t>p{t)))^. 



(5) 



aP 



A frequently used simplification is possible if dissipa- 
tive effects are relevant only on time scales longer than 
any 27r?i/(eQ — ep — e^' -f epi). Then, one can employ a 
full rotating-wave approximation (RWA) for which the di- 
agonal and off-diagonal density matrix elements decouple, 
such that the quantum attractor becomes diagonal, i.e., 

QaP - Oaa^aP- 

Symmetries. — Before evaluating the ratchet current 
we determine two symmetry conditions under which the 
current vanishes. The first one is the generalized parity 
S : {x,t) -^ {-x,t + T 12) [26], which is present if the 
potential and the driving field fulfill the relations 



u{-x)=u[x); E(t + r/2) = -E{t). 



(6) 



Then, the Floquet states obey <f)a{—x,t + T 12) = 
o'a4>a{x,t)^ where CTq = ±1 according to the generalized 
parity. As a consequence, we find pap — —o'aO'pPap and 
thus in particular paa = 0. This means that all Floquet 
states are non-transporting on time-average [17], such that 
any non-vanishing current must stem from off-diagonal 
density matrix elements. The master equation (3) inherits 
a symmetry from the position matrix elements for which 
the relation Xapn = {—^y^^^cra'yf:iXapn holds [12]. This 
leads to the conclusion that the asymptotic state obeys 
g'^p = GaCTpg'^p- Inserting these symmetry relations into 
expression (5) yields J = —J, which implies that the 
ratchet current vanishes in the presence of generalized par- 
ity. 

A second relevant symmetry is time-reversal symmetry 
t -^ —t. It has the consequence that if g{t) is a solution 
of the master equation (3), then g{—t) is a solution as 
well. Time-reversal symmetry is present if the driving field 
obeys 

E{t + Q^E{-t + U), (7) 
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with some appropriate time shift tg, and obviously can 
persist only in the Hamiltonian limit 7 = 0, for which 
the dynamics depends on the initial conditions. In the 
Hamiltonian limit, a meaningful ratchet current requires 
averaging over all possible initial conditions. 

Let us again emphasize that all in the presence of either 
of the two symmetries in eqs. (6), (7), Floquet states are 
non-transporting such that their current paa vanishes, cf. 
Ref. [17]. Since within full RWA by construction g'^p = 
for a ^ (3, while paa = 0, the current (5) vanishes within 
this approximation as well. This in turn means that the 
purely dissipation-induced ratchet current studied below 
can be obtained only from the full master equation (3). 

In order to observe a quantum ratchet current, we need 
to specify the periodic potential u{x) and the driving field 
E{t) such that at least one of the conditions in eq. (6) is 
violated.^ One possibility would be to use a non-reflection- 
symmetric static potential u{x) together with a sinusoidal 
driving [30]. Here, by contrast, we consider a symmetric 
potential and a bichromatic driving field, i.e., 

u{x) = cos(x); E{t) = El cos{ujt)+E2 cos{2Lut + 9), (8) 

which breaks generalized parity provided that both Ei and 
£'2 are non-zero [17]. Indeed, if either _Ei = or _E2 = 0, 
the ratchet current vanishes, both classically and quan- 
tum mechanically. The phase lag 6 allows one to control 
the time-reversal symmetry: If is a multiple of tt, the 
driving field obeys the symmetry condition (7). For any 
other phase lag and Ei,E2 7^ 0, time-reversal symmetry 
is broken, as can be seen in the Husimi functions [31] of 
the Floquet states depicted in fig. 1. Moreover, since the 
transformation (a;, 6) — > {—x, 0±7r) leaves the Hamiltonian 
(2) invariant while it inverts the current, we find that in 
the Hamiltonian limit, the ratchet current obeys [17] 



J{9) = -J{d + TT) =-J{-t 



(9) 



Notably, this relation docs not hold for finite dissipation 
strength 7 > 0. 

Quantum ratchet current and quantum attrac- 
tor. — We start out our numerical studies of the dissi- 
pative quantum system by validating our master equation 
approach for the quantum-classical correspondence. Al- 
ready a classical Brownian particle in the driven periodic 
potential (8) exhibits a rather rich dynamics, ranging from 
regular limit cycles to chaotic motion on strange attrac- 
tors, see fig. 2. The corresponding quantum dynamics is 
even more complex: In the deep quantum regime h > 1, 
it is restricted to a few Floquet states. In the semiclas- 
sical regime ?i <C 1, by contrast, many levels play a role 
and, thus, we expect to find in the Husimi representation 



^For the time evolution of a localized atom cloud in an extended 
periodic optical potential [13,30], one must consider the whole range 
of quasimomenta, k G [—tt,it] [15]. Then the symmetry transforma- 
tions involve eigenstates with opposite k. Neverhteless, the skewed 
initial conditions generally yield a nonzero current even in symmetric 
potentials. 



(X 




Fig. 1: (Color online) Husimi representation of the eigenstate 
for the potential (8) for phase lag — (a) and 8 — — 7r/2 (b). 
The corresponding momenta are paa = (a) and paa ~ —0.15 
(b) (in units of the recoil momentum). The eigenstate has 
been tracked along the corresponding quasienergy band. The 
parameters are h — 0.5, Ei = 1.6, E2 = 2, uj = 1. 



of the density operator signatures of the classical phase- 
space structure. This represents a demanding requirement 
for our master equation formalism. In order to emphasize 
the power of the Floquet master equation (3), we plotted 
the Husimi function of the quantum attractor for both 
a regular limit cycle [fig. 2(a)] and a strange attractor 
[fig. 2(b)]. Comparison with the corresponding classical 
attractors underlines that our formalism is able to cope 
with the semiclassical limit. ^ 

Since the classical attractor shown in fig. 2(a) is 
bounded, it is non-transporting, Jd — 0. The quantum 
attractor, by contrast, supports a very small, but finite dc 
current, Jqm — 0.0025. Note that the dc current for the 
chaotic attractor shown in fig. 2(b) is much larger, namely 
Jci = 0.45 and Jqm — 0.32, respectively. 

Figure 3(a) depicts the ratchet current as a function 
of the phase lag 6 for different dissipation strengths. In 
the Hamiltonian limit 7 = 0, the current"^ vanishes at the 
symmetry points = 0, tt as discussed above. Moreover, 
it complies with relation (9). For finite dissipation, the 
current exhibits multiple current reversals upon chang- 
ing the phase lag 6. This feature is already very pro- 
nounced for 7 = 5- 10"^, which emphasizes that even 
very weak dissipation changes the behavior significantly. 
For much stronger dissipation, 7 — 0.02, the magnitude 
of the current changes slightly, while we still observe sim- 
ilar current reversals. For classical dissipative ratchets, 
such current reversals have been attributed to tangent bi- 
furcations when going from limit cycles towards strange 
attractors [32]. In the Hamiltonian limit of the quantum 
dynamics, however, these bifurcations are absent. There- 
fore, the quantum current reversals may be attributed to 
dissipation as well. 

We next focus on the symmetry point 9 — 0, where the 



^The classical dissipative equations of motion corresponding to 
the Hamiltonian (2) read x = p, p = — 7P + sin(3;) -|- E(t) + 7^1(4). 

^Following Ref. [17], we used in the Hamiltonian limit the initial 
condition ^p{x,to) = (27r)~^", i.e., the zero-momentum plane wave, 
and average the current over the initial time to of the driving. 
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Fig. 2: (Color online) Husimi representation of the quantum 
attractor at stroboscopic times with the corresponding classical 
attractor superimposed (white dots) for E\ = E2 = 2, 'y = 0.1, 
/cbT = 0.1. Driving frequency, phase lag, and effective Planck 
constant are u; = 2, 8 = 7r/2, h = 0.1 (a) and to = 0.87, 
6 = —n/2, h = 0.2 (b). The white line in panel (a) marks the 
corresponding limit cycle. The numerical integration has been 
performed with 45 basis states and e = 0.3. 



Hamiltonian system (2) is time-reversal sym.metric and 
non-transporting. Time-reversal symmetry implies invari- 
ance under p -^ —p, which is perceivable in the Husimi 
representation of the Floquet states at stroboscopic times 
shown in fig. 1(a). Finite dissipation, however, destroys 
time-reversal symmetry, such that the attractor looses the 
symmetry p -^ —p, see fig. 3(b,c), despite the fact that 
it is composed of symmetric, non-transporting Floquet 
states. This reveals that genuine quantum coherence, i.e. 
off-diagonal density matrix elements, play a crucial role for 
both the shape of the attractor and the ratchet current. 
Figure 4(a) shows the dependence of the ratchet current 
on the dissipation strength 7. Both for 7 = and in the 
limit 7 — > 0, the current vanishes. For 7 > we observe a 
purely dissipation-induced quantum ratchet current. This 
current is negative for faint dissipation, but crosses zero 
and becomes positive with increasing dissipation. This 
current reversal behavior resembles the one found for the 
corresponding classical problem [33], but even there has 
not been explained analytically. The dependence of the 
ratchet current on the amplitude of the second harmonic, 
E2 in (8), is shown on Fig. 4(b). 

Approaching the Hamiltonian limit. For dis- 

sipative chaotic quantum systems, the limit of vanishing 
dissipation, 7-^0, deserves some attention. This will 
reveal an intriguing difference between classical and quan- 




Fig. 3: (Color online) (a) Average current J (in units of the re- 
coil momentum) as a function of the phase lag for dissipation 
strength 7 = 0.02 (solid line) and 7 = 5 x 10"* (dashed) com- 
pared to the Hamiltonian case (dash-dotted). The inset shows 
the current for ^ = as a function of the dissipation strength 
7. Lower panel: Husimi representation of quantum attractors 
at stroboscopic times t = nT for 9 — Q: (b) 7 = 5- 10~* and 
(c) 7 = 0.02. The parameters are h — 1, Ei — 1.6, E2 — 2, 
u! = 1, and ksT — 0.5. The numerical integration has been 
performed with 21 basis states and e = 0.3 



tum ratchets. Already in the presence of arbitrarily small 
dissipation, the classical phase-space flow converges to a 
generally unique attractor, while a Hamiltonian system in 
principle preserves memories to its initial state for an ar- 
bitrarily long time. Thus, for 7 = 0, the classical current 
generally depends on the initial conditions, in particular, 
if the initial phase-space distribution overlaps with reg- 
ular manifolds [34] . If the initial condition lies within a 
chaotic manifold, however, the system will in the long- 
time limit be distributed uniformly over the chaotic man- 
ifold and, thus, will eventually become independent of the 
initial condition [11,15]. For the corresponding quantum 
system, this ergodicity is not found owing to the linearity 
of the Schrodinger equation in combination with the finite 
number of states involved. As a consequence, the limit 
of arbitrarily weak dissipation, 7 — > 0, may for quantum 
systems be different from the Hamiltonian case, 7 = 0, as 
we observe in fig. 3(a). 

Conclusions. — We have studied the quantum 
ratchet effect in the weakly dissipative regime in which 
the quantum coherence suffers from decoherence and re- 
laxation. A central property of the corresponding, unique 
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Fig. 4: (Color online) (a) Average current J (in units of the 
recoil momentum) as a function of (a) the dissipation strength 
7 (E-i = 2, ^ = 0) and (b) as a function of the second har- 
monic amplitude, E2 in (8) (7 = 0.02, Q = 7r/2). The other 
parameters are the same as in Fig. 3. 



quantum attractor is a quantum ratchet current, given 
by the time-averaged momentum expectation value. We 
found that even for very weak dissipation, the current dif- 
fers strongly from its corresponding Hamiltonian counter- 
part. The presence of any of two symmetries, namely 
generalized parity and, in the Hamiltonian limit, time- 
reversal symmetry, inhibits a ratchet current. For bichro- 
matic driving, the phase lag Q between the two harmonics 
determines whether the Hamiltonian part is time-reversal 
symmetric or not. If time reversal holds, the current van- 
ishes in the Hamiltonian limit, while for finite dissipation, 
we observe a purely dissipation-induced quantum ratchet 
current which, moreover, possesses current reversals as a 
function of the dissipation strength. 

For cold atoms, the resulting currents are of the order 
10-30% of the recoil momentum, being measurable with 
present experimental techniques [14]. Our study provides 
evidence that cold atoms in driven periodic potentials are 
a natural candidate for studying the complex dynamics 
originating from an intriguing interplay of nonlinearity, 
weak quantum dissipation, and spatio-temporal symmetry 
violation. 
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